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The RENO experiment has observed the disappearance of reactor electron antineutrinos, con- 
sistent with neutrino oscillations, with a significance of 4.9 standard deviations. Antineutri- 
nos from six 2.8 GWth reactors at the Yonggwang Nuclear Power Plant in Korea, are detected 
by two identical detectors located at 294 m and 1383 m, respectively, from the reactor ar- 
ray center. In the 229 day data-taking period between 11 August 2011 and 26 March 2012, 
the far (near) detector observed 17102 (154088) electron antineutrino candidate events with a 
background fraction of 5.5% (2.7%). The ratio of observed to expected numbers of antineutri- 
nos in the far detector is 0.920 ± 0.009(stat.) ± 0.014(syst.). From this deficit, we determine 
sin^ 2^13 = 0.113 ± 0.013(stat.) ± 0.019(syst.) based on a rate-only analysis. 

PACS numbers: 14.60.Pq, 29.40.n, 28.50.Hw, 13.15.+g 



We report a definitive measurement of the neutrino os- 
cillation mixing angle, ^13, based on the disappearance of 
electron antineutrinos emitted from six reactors. Of the 
three mixing angles in the Pontecorvo-Maki-Nakagawa- 
Sakata matrix [1,, , ^13 remains the most poorly known. 
Previous attempts at measuring ^13 via neutrino oscilla- 
tions have obtained only upper limits [3—9]; the CHOOZ 
[3,] and MINOS Q experiments set the most stringent 
limits: sin^ 26*13 < 0.15 (90% C.L.). Recently, indica- 
tions of a non-zero ^13 value have been reported by two 
accelerator appearance experiments, T2K and ML 
NOS [ll|, and by the Double Chooz reactor disappear- 
ance experiment Global analyses of all available 
neutrino oscillation data have indicated central values of 
sin^ 2013 that are between 0.05 and 0.1 (see e.g. 13. l3|). 



During the preparation of this paper, the Daya Bay ex- 
periment reported the measurement of a non-zero value 
for 6*13 M- 

Reactor experiments with a baseline distance of ^1 km 
can neglect the disappearance of Pe driven by 612 and 
Am|i and, thus, unambiguously determine the mixing 
angle 6*13 based on the survival probability of electron 
antineutrinos, 



survival 



1 - sin^ 26ii3 sin2(1.267ATO^iL/£;), (1) 



where E is the energy of antineutrinos in MeV, and 
L is the baseline distance in meters between the reac- 
tor and detector. The well measured value of Am|2 = 
{2.32toll) X 10-3 eV2 can be substituted for Am^i 
in Eq. (1). 

The detection methods and setup of the RENO ex- 
periment are discussed in detail elsewhere 17[. In this 
Letter, only the components relevant to this measure- 
ment are reviewed. Two identical antineutrino detectors 
are located at 294 m and 1383 m, respectively, from the 
center of reactor array to allow a relative measurement 
from a comparison of the measured neutrino rates. The 
measured far-to-near ratio of antineutrino fluxes can con- 
siderably reduce systematic errors coming from uncer- 
tainties in the reactor neutrino flux, target mass, and 
detection efficiency. The relative measurement is inde- 
pendent of correlated uncertainties and helps minimize 
uncorrelated reactor uncertainties. The far (near) detec- 
tor is under a 450 (120) meters of water equivalent rock 
overburden. 

Six pressurized water reactors, each with maximum 
thermal output of 2.8 GWth (reactors 3, 4, 5, and 6) 
or 2.66 GWth (reactors 1 and 2), are situated in a line 
with roughly equal spacings and span a total distance of 
^1.3 km. The positions of the two detectors and the six 
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reactors were surveyed with GPS and total station to de- 
termine the basehne distances between the detectors and 
reactors to accuracies better than 10 cm. The reactor- 
flux weighted basehne is 408.56 m for the near detector, 
and 1443.99 m for the far detector. 

The reactor is detected via the inverse beta decay 
(IBD) reaction, S'e + P e+ + n. Detectors based on 
hydrocarbon hquid scintiUator (LS) provide free protons 
as a target. The coincidence of a prompt positron signal 
and a delayed signal from neutron capture by Gadolinium 
(Gd) provides the distinctive IBD signature. 




FIG. 1. A schematic view of a RENO detector. The near and 
far detectors are identical. 

Each RENO detector (Fig. [ij consists of a main inner 
detector (ID) and an outer veto detector (OD). The main 
detector is contained in a cylindrical stainless steel vessel 
that houses two nested cylindrical acrylic vessels [ll|. The 
innermost acrylic vessel holds the 18.6 m'^ (16 t) ~0.1% 
Gd-doped LS as a neutrino target. It is surrounded by a 
7-catcher region with a 60 cm thick layer of Gd-unloaded 
LS inside an outer acrylic vessel. Outside the 7-catcher 
is a 70 cm thick buffer region filled with 65 tons of min- 
eral oil. Light signals emitted from particles interacting 
in ID are detected by 354 10-inch Hamamatsu R7081 
photomultiplier tubes (PMTs) that are mounted on the 
inner wall of the stainless steel container. The 1.5 m 
thick region of the OD that is external to the ID is filled 
with highly purified water. The OD is equipped with 67 
10-inch R7081 PMTs mounted on the wall of the veto 
vessel. The LS is developed and produced as a mixture 
of linear alkyl benzene (LAB), PPO, and bis-MSB. A 
Gd-carboxylate complex using TMHA was developed for 
the best Gd loading efficiency into LS and its long term 
stability [19 1. 



Event triggers are formed by the number of PMTs 
with signals above a '^Q.S photoelectron (p.e.) thresh- 
old (NHIT). An event is triggered and recorded for an 
IBD candidate if the ID NHIT is larger than 90, corre- 
sponding to 0.5~0.6 MeV and weU below the 1.02 MeV 
minimum energy of an IBD positron signal. The ID trig- 
ger provides no loss of IBD candidates. 

The detectors are calibrated using radioactive sources 
and cosmic-ray induced background event samples. Ra- 
dioisotopes of ^^''Cs, ^^Ge, ^°Co, and ^^^Cf are periodi- 
cally deployed in the target and 7-catcher by a step mo- 
torized pulley system in a glove box. The detectors' en- 
ergy response stability is continuously monitored using 
cosmic-ray produced neutron captures on H and Gd. 

The event energy is determined from the total p.e.- 
charge {Qtot) that is collected by the PMTs, corrected 
for gain variation. The energy calibration constant of 
250 p.e. per MeV is determined from the peak ener- 
gies of various radioactive sources deployed at the cen- 
ter of the target. The obtained energy resolution is 
(5.9/^£;(MeV) + 1.1)%, common for both detectors. 

In this analysis, an IBD event requires a delayed signal 
from a neutron capture on Gd and, thus, the fiducial 
volume naturally becomes the entire target vessel region 
without any vertex position cuts. There is some spill-in 
of IBD events that occur outside the target and produce 
a neutron capture on Gd in the target, which enhances 
the detection efficiency. 

The following criteria are applied to select IBD can- 
didate events: (1) Qmax/Qtot < 0.03 where Qmax is the 
maximum charge of a PMT, to eliminate PMT flasher 
events and external 7-ray events; (2) a cut rejecting 
events that occur within a 1 ms window following a cos- 
mic muon traversing the ID with an energy deposit {E^j) 
that is larger than 70 MeV, or with iJ^ between 20 MeV 
and 70 MeV for OD NHIT > 50; (3) events are rejected 
if they are within a 10 ms window following a cosmic 
muon traversing the ID if is larger than 1.5 GeV; (4) 
0.7 MeV < Ep < 12.0 MeV; (5) 6.0 MeV < Ed < 12.0 
MeV where Ep {Ed) is the energy of the prompt (de- 
layed) event; (6) 2 /zs < Afe+„ < 100 jis where Aig+„ is 
the time difference between the prompt and delayed sig- 
nals; (7) a multiplicity requirement rejecting correlated 
coincidence pairs if they are accompanied by any preced- 
ing ID or OD trigger within a 100 /iS window before their 
prompt candidate. 

Applying the IBD selection criteria yields 17102 
(154088) candidate events or 77.02±0.59 (800.8±2.0) 
events/day for a live time of 222.06 (192.42) days in 
the far (near) detector. In the final data samples, some 
uncorrelated (accidentals) and correlated (fast neutrons 
from outside of ID, stopping muon followers, and /3-n 
emitters from ^Li/*He) background events survive the 
selection requirements. 

The uncorrelated background is due to accidental co- 
incidences from the random association of a prompt-like 
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TABLE I. Event rates of the observed candidates and the 
estimated background. 



Detector 


Near 


r ar 


Selected events 


154088 


17102 


Total background rate (per day) 


21.75±5.93 


4.24±0.75 


IBD rate after background 


779.05±6.26 


72.78±0.95 


subtraction (per day) 






DAQ Live time (days) 


192.42 


222.06 


Detection efficiency (e) 


0.647±0.014 


0.745±0.014 


Accidental rate (per day) 


4.30±0.06 


0.68±0.03 


®Li/*He rate (per day) 


12.45±5.93 


2.59±0.75 


Fast neutron rate (per day) 


5.00±0.13 


0.97±0.06 



event due to radioactivity and a delayed-like neutron cap- 
ture. The remaining rate in the final sample is estimated 
by measuring the rates of prompt- and delayed-like events 
after applying all the selection criteria other than (6), and 
calculating the probability of random association in the 
At window for IBD selection, leading to 4.30±0.06 (near) 
or 0.68±0.03 (far) events per day. 

The ^Li/®He /3-n emitters are mostly produced by en- 
ergetic muons because their production cross sections in 
carbon increase with muon energy [20—22]. The back- 
ground rate is estimated from a sample prepared by a de- 
layed coincidence between an energetic (i?^ > 0.5 GeV) 
muon and the following IBD-like pair of events. The 
^Li/®He /3-n background rate in the final sample is ob- 
tained as 12.45±5.93 (near) or 2.59±0.75 (far) events per 
day from a fit to the delay time distribution with an ob- 
served mean decay time of ~250 ms. 

An energetic neutron entering the ID can interact in 
the target to produce a recoil proton before being cap- 
tured on Gd. Fast neutrons are produced by cosmic 
muons traversing the surrounding rock and the detector. 
The background rate is estimated by extrapolating the 
energy spectral shape of events with 12 MeV < i5p < 30 
MeV, to the IBD signal region, assuming a flat spectrum 
of the fast neutron background. The estimated fast neu- 
tron background is 5.00±0.13 (near) or 0.97±0.06 (far) 
events per day. The total background rate is estimated 
to be 21.75±5.93 (near) or 4.24±0.75 (far) events per day 
and summarized in Table ID 

Both the prompt energy and flasher requirements are 
almost fully (99.8%) efficient. The fraction of neutron 
captures on Gd is evaluated to be (85.5±0.7)% using MC 
and Cf source data. The Ate+„ requirement efficiency 
is determined to be (92.1±0.5)% from MC and data. The 
fraction of neutron captures on Gd that satisfy the 6.0 
MeV threshold requirement is (95.2±0.5)%. The overah 
efficiency for finding a delayed signal as an IBD candidate 
pair is (74.9±1.0)%. The spill- in IBD events result in a 
net increase in the detection efficiency of 2.2%. The com- 
mon detection efficiency is estimated to be (76.5±1.4)% 
using a Monte Carlo simulation (MC) and data. 



The fractional losses of IBD events due to the muon 
veto are determined to be (11.30±0.04)% (near) or 
(1.36±0.02)% (far), by summing the time spent in ve- 
toing events after muons. The fractional losses of IBD 
events due to the multiplicity cut is calculated to be 
(4.61±0.04)% (near) or (1.22±0.07)% (far), based on the 
ID trigger rate and the veto window from an IBD prompt 
candidate. The efficiencies for detecting IBD events are 
found to be (64.7±1.4)% (near) and (74.5±1.4)% (far). 

The absolute uncertainties of the efficiencies are corre- 
lated between the two detectors. Only differences be- 
tween the two identical detectors are taken as uncor- 
related uncertainties. The systematic uncertainties are 
summarized in Table II. 

TABLE II. Systematic uncertainties in the reactor neutrino 
detection. 



Reactor 





Uncorrelated 


Correlated 


Thermal power 


0.5% 




Fission fraction 


0.7% 




Fission reaction cross section 




1.9% 


Reference energy spectra 




0.5% 


Energy per fission 




0.2% 


Combined 


0.9% 


2.0% 


Detection 




Uncorrelated 


Correlated 


IBD cross section 




0.2% 


Target protons 


0.1% 


0.5% 


Prompt energy cut 


0.01% 


0.1% 


Flasher cut 


0.01% 


0.1% 


Gd capture ratio 


0.1% 


0.7% 


Delayed energy cut 


0.1% 


0.5% 


Time coincidence cut 


0.01% 


0.5% 


Spill-in 


0.03% 


1.0% 


Muon veto cut 


0.02% 


0.02% 


Multiplicity cut 


0.04% 


0.06% 


Combined (total) 


0.2% 


1.5% 



Uncorrelated relative uncertainties are estimated by 
comparing the two detectors. The IBD differential cross 
section is taken from Ref. 23]. The total number of 
free protons in the target is 1.189 x 10^° with an uncer- 
tainty of 0.5%, determined from measurements of the LS 
weight and composition. The relative energy scale dif- 
ference between the detectors is determined to be 0.2% 
from comparison of the peak energy values for several 
radioactive calibration sources, IBD delayed events, and 
cosmic muon induced spallation-neutron captures on H 
and Gd. The energy scale difference is found to corre- 
spond to a relative uncertainty in the efficiency of the 
delayed energy of 0.1% using data. The Gd-LS was com- 
monly produced and then divided equally and filled into 
the two detectors to ensure that the Gd concentration 
and the target protons of the near and far detectors are 
identical. This procedure for filling the targets results 
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in a difference in the number of the target protons that 
is less than 0.1%. The difference in the measured neu- 
tron capture time between the detectors is less than 0.2 
IIS, corresponding to Gd concentration differences of less 
than 0.1%. The relative uncertainty of Gd capture ra- 
tio is less than 0.1% accordingly. The remaining relative 
uncertainties are close to 0.01%, and the combined un- 
certainty common to the both detectors is 0.2%. A more 
detailed discussion on the systematic uncertainties will 
be presented in a future publication. 

The antineutrino flux depends on thermal power, fis- 
sion fractions of the four isotopes, energy released per 
fission, and fission and capture cross-sections. The un- 
certainty associated with the thermal power, provided 
by the power plant, is 0.5% per core and fully correlated 
among the reactors 2J]. The relative fission contribu- 



tions of the four main isotopes are evaluated for the fuel 
cycle with 4^10% uncertainties, using the Westinghouse 
ANC reactor simulation code [25|. The uncertainties of 
the fission fraction simulation contribute 0.7% of the j>e 
yield per core to the uncorrelated uncertainty. The as- 
sociated antineutrino fiux is computed based on the De 
yield per fission [26|] and the fission spectra [27—31], lead- 
ing to a 1.9% correlated uncertainty that has little effect 
on the 9i3 determination. The thermal energy released 
per fission is given in Ref. [s^ , and its uncertainty results 
in a 0.2% correlated uncertainty. We assume a negligible 
contribution of the spent fuel to the uncorrelated uncer- 
tainty in this analysis. 

All reactors were mostly in steady operation at the full 
power during the data-taking period, except for reactor 2 
(R2), which was off for the month of September 2011, and 
reactor 1 (Rl), which was off from February 23 2012 for 
fuel replacement. Figure [2] presents the measured daily 
rates of IBD candidates after background subtraction in 
the near and far detectors. The expected rates assuming 
no oscillation, obtained from the weighted fiuxes by the 
thermal power and the fission fractions of each reactor 
and its baseline to each detector, are shown for compar- 
ison. 

The ratio of measured to expected events in the far 
detector is 

R = 0.920 ± 0.009(stat.) ± 0.014(syst.), 

which indicates a clear deficit. To determine the value 
of sin^ 26'i3 from the deficit, a with pull terms on the 
correlated systematic uncertainties 



is used. 



E 



[Ks + 6d - (1 + « + Crf) + fr)Niil 



d=N,F 



obs 



6 

E 



(2) 




where d is an index denoting the near detector (N) or 
the far detector (F), r corresponds to reactors 1 through 



Aug 29 Oct 28 Dec 27 Feb 25 
2011 „ , 2012 
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FIG. 2. Measured daily-average rates of reactor neutrinos 
after background subtraction in the near and far detectors as 
a function of running time. The solid curves are the predicted 
rates for no oscillation. 



6, A^obs is the number of observed IBD candidates in 
each detector after background subtraction, and N^^p is 
the number of expected neutrino events, including de- 
tection efficiency, neutrino oscillations, and contribution 
from the r-th reactor to each detector determined from 
baseline distances and reactor fluxes. A global normal- 
ization a is taken free and determined from the fit to the 
data. Then, a is constrained by the normalization uncer- 
tainty of 2.5%, coming from correlated uncertainties, to 
the value obtained from the fit. The uncorrelated reac- 
tor uncertainty is 0.9% (iTr), the uncorrelated detection 
uncertainty is 0.2% (u^), as listed in Table HIl and cr^ is 
the background uncertainty listed in Table HI /r, £,d, and 
bd are corresponding pull parameters. 
The best-fit value thus obtained is 

sin^ 2013 = 0.113 ± 0.013(stat.) ± 0.019(syst.), (3) 

and excludes the no-oscillation hypothesis at the 4.9 stan- 
dard deviation level. 

Figure [3] shows the distribution as a function of 
sin^ 2013 1 and the ratios of the measured reactor neutrino 
events, relative to what is expected without oscillation at 
both detectors. We observe a clear deficit of 8.0% for the 
far detector, and of 1.2% for the near detector, conclud- 
ing a definitive observation of reactor antineutrino disap- 
pearance consistent with neutrino oscillations. The sur- 
vival probability due to neutrino oscillation at the best-fit 
value is given by the curve. 

The observed spectrum of IBD prompt signals in the 
far detector is compared to non-oscillation expectations 
based on measurements in the near detector in Fig. 2] 
The spectra of prompt signals are obtained after sub- 
tracting backgrounds shown in the inset. The disagree- 
ment of the spectra provides further evidence of neutrino 
oscillation. 

In summary, RENO has observed reactor antineutri- 
nos using two identical detectors each with 16 tons of 
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Weighted Baseline [m] 

FIG. 3. The distribution as a function of sin^20i3. Bot- 
tom: Ratio of the measured reactor neutrino events relative 
to the expected with no oscillation. The curve represents the 
oscillation survival probability at the best fit, as a function of 
the flux-weighted baselines. 



Gd-loaded liquid scintillator, and a 229 day exposure to 
six reactors with total thermal energy 16.5 GWt^- In the 
far detector, a clear deficit of 8.0% is found by compar- 
ing a total of 17102 observed events with an expectation 
based on the near detector measTircmcnt assTiming no os- 
cillation. From this deficit, a rate-only analysis obtains 
sin^ 2013 = 0.113 ± 0.013(stat.) ± 0.019(syst.). The neu- 
trino mixing angle Q^?, is measured with a significance of 
4.9 standard deviation. 
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